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ABSTRACT 

We report the parallax and proper motion of millisecond pulsar J0030+0451, one of thirteen known isolated 
millisecond pulsars in the disk of the Galaxy. We obtained more than 6 years of monthly data from the 305 
m Arecibo telescope at 430 MHz and 1410 MHz. We measure the parallax of PSR J0030+0451 to be 3.3 ± 
0.9 mas, corresponding to a distance of 300 ± 90 pc. The Cordes and Lazio (2002) model of galactic electron 
distribution yields a dispersion measure derived distance of 317 pc which agrees with our measurement. We 
place the pulsar's transverse space velocity in the range of 8 to 17 km s" 1 , making this pulsar one of the slowest 
known. We perform a brief census of velocities of isolated versus binary millisecond pulsars. We find the 
velocities of the two populations are indistinguishable. However, the scale height of the binary population is 
twice that of the isolated population and the luminosity functions of the two populations are different. We 
suggest that the scale height difference may be an artifact of the luminosity difference. 

Subject headings: binaries — pulsars: individual(J0030+045 1 ) — solar neighborhood — solar wind — stars: 
distances 



1. INTRODUCTION 

A pulsar parallax can be combined with a measurement 
of the pulsar's dispersive delay (the Dispersion Measure or 
DM) to provide an accurate measure of the free electron 
density along the line of sight (LOS). PSR J0030+0451 
is one of fewer than a doz en pulsars to have its paral- 
lax measured via timing dKasp i. Tavlor. & Rvbal [19941 
Camilo. Foster & Wolszczanl Il994[ ISandhu et all Il997t 
Toscano et all Il999bt IWolszczan et al. I2000at iTacobv et al l 
20031 iHo tan. Baile s. & Ordl 120041 lLohmer et alJ 2004; 
Splave ret alJ |2005). Another dozen have been measured via 
VLBI (e.g. see Brisken et al. 2002, Chatterjee et al. 2004, and 
also this URL 1 ). These measurements are important because 
they give us most of o ur knowledge about the galactic thermal 
electro n distribution dCord es & La ziol 120021 iToscano et all 
19995 lTavlo7 Manches ter. & Lvnell993l) . 

In addition, PSR J0030+045 1 presents a rare evolutionary 
case as an isolated millisecond pulsar (MSP). In the most pop- 
ular MSP evolutionary model, MSPs are formed via accretion 
of matter from a companion star. The incoming matter adds 

1 http://www.astro.cornell.edu/~shami/psrvlb/parallax.html 



to the pulsar's angular momentum, i.e., the pulsar is "spun 
up." However, isolated MSPs present a conundrum: they were 
presumably spun up, yet they are without a companion which 
would have done so. One possibl e scenario is that the pul- 
sar ha s ablated its companion ( Ruderman, Shaha m7& Tavanil 
1989). 

We expect MSPs to have lower velocities than the reg- 
ular population, because the kick from the supernova pro- 
genitor had to be small enough to leave the binary intact. 
To produce an isolated MSP, the binary must remain intact 
during and after the supernova, but then after the spin-up 
phase the companion must leave the system or be evaporated. 
Several authors have debated whether isolated MSP veloci- 
ties are lower, higher, or indistinguishable from those of the 
general population of MSPs. McLaughlin et al. (2004) sug- 
gest we might expect isolated MSPs to have higher veloci- 
ties. They argue that if isolated MSPs are formed by abla- 
tion, we would expect them to form from the tighter bina- 
ries which are more susceptible to ablation. The correlation 
between tight b i naries and higher velocities is su ggested by 
iTauris & Bailesl (119961) . iMcLaughl in et alJ (|2004) present the 
argument for faster velocities for isolated MSPs as a counter- 
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point to their timing proper motion and scintillation measure- 
m ents which suggest the opposite, as do the measurements 
oflJohnston. Nicastro 1 & Koribalski ( 119981) and lToscano et alJ 
( 1999b). Hobbs et alJ fe005l) . however, find the velocities of 
the populations to be indistinguishable. We present a mea- 
surement of the transverse velocity of PSR J0030+045 1 which 
is unusually small, even compared to the isolated MSP popu- 
lation. We reconsider the question of the velocity of isolated 
MSPs as compared to the binary MSP population. 

Owing to its small timing residual, ~1 /is, PSR 
J0030+045 1 is a good candidate for membership in the Pulsar 
Timing Array (PTA), which is a collection of pulsars that will 
be used for detecting g ravitational radiation (Jaffe & Backer 
120031 Uenet et alJ l2004). For this reason, continued refine- 
ment of the timing model of PSR J0030+0451 is impor- 
tant. In fact, the PTA, as it is conceived, is an interferom- 
eter, so a variety of baselines will be important to its oper- 
ation. PSR J0030+0451 may be particularly useful in this 
regard because it has large angular separation from PSRs 
B 1855+09, J 17 13+07 and J043 7-4715, which are among the 
most stable and precise pulsars (Kaspi, Tavlor, & Rvba 1994; 
Ivan Straten et all2001tlLommenl20 01). 

In §2 we present a significant refineme nt to the timing 
model previously published ( Lommen et al. 2000). We dis- 
cuss our measurements of parallax and proper motion in |3] 
We include the effects of the solar wind in our analysis, which 
we discuss in §4. In §5, §6, and §7 we discuss the space veloc- 
ity of J0030+0451, corrections to its measured period deriva- 
tive, and the implications of its measured distance for the local 
interstellar medium (LISM). In §8 we summarize our conclu- 
sions. 

2. ARECIBO OBSERVATION AND DATA REDUCTION 

We have conducted timing observations over a 6.5 year pe- 
riod, from 1997 December to 2004 July, at the Arecibo Obser- 
vatory, using the Arecibo-Berkeley Pulsar Processor (ABPP) 
and the Princeton Mark IV s ystem. Details of AB PP sig- 
nal processing can be found in Lomm en et al. (2000), where 
the first three years of these data are pr esented. Details o f 
Mark IV signal processing can be found in Stairs et al. ( 2000). 

Over the course of each observation, incoming signals from 
two orthogonal polarizations were coherently dedispersed, 
squared to obtain power measurements, and folded modulo 
the pulse period for intervals of three minutes. Opposite po- 
larizations were summed, using absolute flux calibrations of 
the receivers whenever possible or, when absolute flux cali- 
bration was not available, using the system temperature of the 
pulsed noise source as published by telescope operations. 

Pulse times of arrival (TOAs) were derived from the cal- 
ibrated three minute integrations using conventional algo- 
rithms. The TOAs from a given frequency and instrument on 
a given day were averaged into a single effective TOA for that 
day, the uncertainty of which was estimated from the spread 
of the individual three-minute TOAs. 

We performed a weighted fit to the averaged TOAs using 
TEMPO 2 . We used ecliptic coordinates to minimize covari- 
ance between the two components of the position and proper 
motion measurements. This is usually necessary when the 
ecliptic latitude is low (1.44° in our case). Table 1 shows 
updated spin, astrometric, and other parameters for PSR 
J0030+0451. The root mean squared (rms) of the timing 
residuals quoted in Table 1, is calculated using TOAs from 

2 See http://pulsai-.princeton.edu/tempo 




51000 52000 53000 
Day 



FIG. 1 . — Residual vs day for the model presented in Table HI Circles are 
ABPP data. Squares are Mark IV data. Filled points are 1410 MHz. Open 
points are 430 MHz data. 



profiles averaged over 30 minutes. We estimated uncertain- 
ties in all the fitted parameters by doubling errors given by 
TEMPO, an ad hoc procedure that attempts to account for tim- 
ing noise and other possible systematic errors. The quoted 
uncertainty in parallax also includes a term added in quadra- 
ture (0.3 mas) due to the solar wind (see |j4}. Residual arrival 
times after subtracting the best fit are shown in Figure^ 

We examined timing solutions in which DM was allowed 
to vary over time, and we found that they did not significantly 
improve the fit quality or change the timing model parameters. 
However, it is interesting to compare our upper limit on DM, 
the time derivative of DM , to the cor relation between DM and 
DM measured by Ba cker et aT](ll993l) . When DM is included 
in the timing model, its best-fit value is (1.9 ± 1.8) x 10~ 5 
cm~ 3 pc yr" 1 , which implies an upper limit (with 95% con- 
fidence) of 5.5 x 10~ 5 cm~ 3 pc yr" 1 . This va lue is smaller than 
expected according to Backer et al. ( 1993) by about a factor 
of 5. If more small-DM small-DM pulsars are found it might 
suggest that DM is a stronger function of DM than we would 
expect for wedgelike thermal plasma perturbations distributed 
randomly along the LOS. 

3. PROPER MOTION AND PARALLAX 

Evidence for a significant parallax measurement is shown in 
Figure |2] which compares residuals of the pulse arrival times 
with and without parallax incorporated into the timing model. 
The figure shows averaged timing residuals versus day num- 
ber binned in increments of 18 days and folded over a half- 
year period to provide the best visibility for the parallax sig- 
nature. The data have been fit for proper motion but not for 
parallax in the left part of Figure [2] Data on the right have 
been fit for both proper motion and for parallax. We have 
superimposed the best-fit parallax curve onto the pre-fit data. 
On close examination of Figure |2] there is a subtle difference 
in the positions of the data points relative to the fit curve (left 
plot) and the horizontal axis (right plot). This is not an unex- 
pected result of (a) a global fit to all the parameters with and 
without parallax and (b) the post-fit binning. 

Our best-fit value for parallax, 7r = 3.3 ± 0.9 mas, includes 
both measurement uncertainty of 0.88 mas and systematic un- 
certainty due to the solar wind model, 0.3 mas, which will be 
discussed in the following section. 

We measure a proper motion of -5 .74 ±0.09 mas yr" 1 in the 
plane of the ecliptic. Proper motion out of the ecliptic plane 
is naturally difficult to measure given the pulsar's position, a 
problem which was compounded by its high covariance with 
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TABLE 1 

Parameters for pulsar J0030+045 1 



Parameter Value" 



Ecliptic longitude (deg) 8.91036695(6) 

Ecliptic latitude (deg) 1.445692(3) 

Period (s) 0.00486545320829334(3) 

Period derivative (ss" 1 ) 1.0162(1) Xl(T 20 

Proper motion in ecliptic longitude (mas yr -1 ) —5.74(9) 

Absolute value of proper motion in ecliptic latitude (mas yr~' ) . . . < 10 

Dispersion measure (pccirT 3 ) 4.3328(2) 

Epoch (MJD) 52035 

Solar n a (e - cm" 3 ) 6.9(2.1) 

Solar ri (e~cnr 3 yr _1 ) 0.0(3) 

Number of epochs of data b 56 

Timing data span (MJD) 50,790 - 53,280 

Right ascension(J2000) 00 30 27.4308(6) 

Declination(J2000) +04 51 39.72(1) 

Galactic longitude (deg) 113.141135(1) 

Galactic latitude (deg) -57.611237(3) 

DM derived distance (pc) c 317(25) 

Characteristic age (yr) 7.8 X 10 9 

Magnetic field strength (G) d 2.7 X 10 8 

Column electron density along LOS (e - cm~ 3 ) 0.014(2) 

Parallax (mas) 3.3(9) 

Parallax derived distance (pc) 300(90) 

Magnitude of transverse velocity (km s ) 8—17 

RMS residual at 430 MHz (ps, ABPP/MarklV) 2.3/1.0 

RMS residual at 1410 MHz (jus, ABPP/MarklV) 2.2/2.2 



"Uncertainties in parentheses refer to the last digit quoted. Note that a iiq of 6.9 cm -3 was 
assumed in all fits. 
b We define "epoch" as data separated by 2 weeks or more. 
c ModelfromNE2001. 
d B = 3.2xl0 19 GPW'o] 1 -' 2 . 
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FIG. 2. — Timing residuals for J0030+0451 folded over a period of a half year and binned every 18 days. Left: Residuals with no parallax fit. Right: Residuals 
after removal of parallax of 3.3 mas from the data in the left hand panel. The curve that represents the 3.3 mas parallax is shown in the left hand panel. Day is 
when the Earth-Sun-pulsar angle is 90° . 



variations in dispersion measure. We quote a conservative up- 
per limit on the magnitude of the out-of -plane proper motion 
of 10 mas yr" 1 . 

We can combine our distance measurement with our proper 
motion measurements to yield the pulsar's velocity. We obtain 
a velocity of -8. 3 ±0.1 km s -1 in the plane of the ecliptic and 
an upper limit of 14.4 km s" 1 out of the plane of the ecliptic. 
Thus, the pulsar's transverse velocity lies bet ween 8 and 17 
km s" 1 . This confirms the value presented by iNicastro et alJ 
fcOOll) . 9±6km s" 1 , which they found using scintillation 
measurements. 

The measured proper motion of PSR J0030+0451 results 
not just from its motion relative to its local standard of rest 



(LSR), but also from the difference between its LSR and the 
LSR of the sun and from the solar motion. After considering 
these effects (see §5 for details), we calculate its transverse 
velocity relative to its LSR to be between 4 and 20 km s" 1 . 

This is one of the slowest transverse velocities measured 
for any pulsar. Young pulsars (those which ha ve not been 
spun up) have a mean velocity of 265 km s" 1 dHobbs et alJ 
2005), a factor of three times faster than MSPs (§5; see also 
Hobbs et al. 2005; Cordes & Chernoff 1997; Nice & Taylor 
1995). PSR J0030+0451 is, in fact, only one tenth as fast as 
the average MSP (§5). 

4. THE SOLAR SYSTEM ELECTRON DENSITY 
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PSR J0030+0451 is roughly in the ecliptic plane (eclip- 
tic latitude is 1.44°) where the solar system provides addi- 
tional disp ersion caused b y the charged particles of the solar 
wind. Ilssautier et al.l ( 120011) show that the solar electron den- 
sity, n e , can be modeled as «o/r 2 , where r is the radial dis- 
tance to the sun in astronomical units, and m is roughly 10 
cm -3 ( Snla ver et alJl2005l Ilssautier et al.ll2.001h . This simple 
model is the default in TEMPO. Maintaining the 1 / r 2 depen- 
dence, we searched x 2 space for the best value of «o, which 
we found to be 6.9 ±2.1 cm" 3 . 

Reducing «o from 9.0 to 6.9 cm" 3 reduces the fitted paral- 
lax by about 0.3 mas which increases the distance by about 30 
pc. These two parameters are covariant because the model of 
the solar wind yields an annual pattern of pulse delays with a 
strong cusp when the pulsar is behind the sun with respect to 
the earth. A Fourier decomposition of this delay pattern yields 
significant terms with annual periodicity (covariant with pul- 
sar position), semi-annual periodicity (covariant with paral- 
lax), as well as higher order terms. The additional 0.3 mas 
uncertainty in parallax due to the solar wind was added in 
quadrature to the doubled TEMPO error. 

Note that in the case of PSR J 17 13+0747, another ecliptic 
plane MSP, ISolaver et alJ (120051) found it necessary to elimi- 
nate data within 30° of the sun. In our case we did not find 
that this significantly improved the fit, nor did it substantially 
change any fitted parameters. 

We wondered if it would be possible to measure a change in 
the solar wind on its 1 1 -year cycle. Using TEMPO we mapped 
X 2 space in the range of -10 cm" 3 yr" 1 < ri„ < 10 cm" 3 yr" 1 
but found that the minimum in \ 2 occurred at an essentially 
null value of 0.0 ± 0.3 cm" 3 yr" 1 , indicating that the change is 
currently beyond our measurement capability. 

5. VELOCITIES OF MSPS: ISOLATED VS BINARY 

Hobbs et al.l (120051) presents an extensive study of the ve- 
locities of various sub-groups of pulsars, including isolated 
and binary recycled pulsars. They find that isolated recy- 
cled pulsars are not significantly slower (77 ± 16km s" 1 ) than 
binary recycled pulsars (89 ± 15 km s" 1 ). When the veloc- 
ity of PSR J0030+0451 (using 9 km s"l) is added to the 
sample the average velocity of isolated recycled pulsars be- 
comes 68 ± 16km s , which is still not significantly dif- 
ferent from th e average of the binary recycled pulsars . In 
contrast, both Johnston. Nicastro. & Koribalski (1998) and 
iToscano et alJ (1999b) previously claimed to see evidence that 
isolated MSPs are slower than binary MSPs. However, when 
one uses a more recent Galactic ele ctron density model to 
estimate the distances to the pulsars (Cordes & Lazio 2002), 
the discrepancy disappears. This is a potent reminder that a 
distance model has a significant e ffect on conclusions drawn 
from it. iMcLaughlin et alJ (120041) also find the isolated MSP 
population to have a slightly lower average velocity (70 
km s" 1 ) as compared to the binary population. They cau- 
tion against drawing conclusions from velocity data which are 
based on imprecise distances. 

We performed an analysis of the velocity data that is 
slightly different from the analysis done by iHobbs et al.l 
(2005) with similar results. We defined an MSP to be any 
pulsar with period P < 0.01 s, which prov ides a narrower sam- 
ple than that used bv lHobbs et alJ ( 120051) who defined MSPs 
to be any pulsar with period P < 0.1 s and period derivative 
P < 10~ 17 . Our criterion means that the binary pulsars in 
our sample nearly all have helium white dwarf companions, 
whereas Hobbs et al. used binaries with a mix of companion 
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FIG. 3. — Histograms of height above the galactic plane for the isolated 
MSP population (upper) and the binary MSP population (lower). 



types. 

The 29 MSPs in the Galactic Disk with measured proper 
motions are shown in Table 13 PSR J 1730-2304, which has 
no measured declination proper motion, has been included in 
the table for completeness but has not been included in any 
of the following calculations. We corrected each pulsar's ve- 
locity to its LSR as follows. We used the measured proper 
motion and distance to calculate a three dimensional vector 
representing the (two dimensional) transverse motion of the 
pulsar in the reference frame of the Sun. We then removed the 
solar motion and rotated the resulting vector from the LSR of 
the sun to the LSR of the pulsar. Finally, we recovered those 
two components of the vector which are perpendicular to the 
line of sight. This computation required selecting a value for 
the unknown LOS velocity of the pulsar; we chose a value 
appropriate for a star at rest in the pulsar's LSR. 

This corrected velocity is listed in the second to last col- 
umn of Table 3. The average corrected velocity of the iso- 
lated MSPs is 86 ± 19 km s" 1 whereas the average corrected 
velocity of all the binary MSPs is 91 ±28 km s" 1 . (The un- 
corrected averages are 79 and 90 km s" 1 respectively). If one 
allows the sample to include only those proper motions which 
have been measured to better than 2a the average corrected 
velocities are 86 ± 19 km s" 1 and 99 ± 33 km s" 1 respectively. 
(Uncorrected averages are 79 and 99 km s" 1 .) In each case 
the isolated MSP population is indistinguishable from the bi- 
nary MSP population. The 2cr cutoff in velocity introduces a 
selection of higher velocity pulsars. Thus, the average veloc- 
ities are higher in that case. 

An alternative statistic for evaluating the dynamics of pul- 
sar populations is the distribution of heights above or below 
the galactic plane, z. For the pulsars listed, one finds that the 
standard deviation from zero for the binary MSP population 
is twice that of the isolated MSP population: 570±90pc vs 
280 ± 65 pc. Figure[3]shows a histogram of z for each popula- 
tion. The isolated MSP population is represented in the upper 
half of the figure, the binary MSP population in the lower half. 

Figure|3]shows that the known isolated MSPs are closer to 
the Plane than are the known binary MSPs. This could be ei- 
ther a reflection of differences in the intrinsic spatial distribu- 
tions of the two types of MSPs, or a selection effect. A smaller 
intrinsic spread in scale heights for isolated MSPs is only pos- 
sible if that population also has a smaller intrinsic velocity dis- 
tribution, so that the objects do not travel as far from the Plane 
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TABLE 2 

Velocities of Millisecond Pulsars in the Galactic Disk 



Pulsar fi a fig Distance v< Reference 

(masyr -1 ) (masyr~') (pc) (kms -1 ) 



Isolated MSPs 



J0O3O+0451 


^ =-5.84 


± 


0.09 


M < io 




310P 


<20 


This work 


J07 11-6830 


-15.7 


± 


0.5 


15.3 


± 0.6 


860" 


113 


Toscano et al. (1999a) 


J1024-0719 


-41 


± 


2 


-70 


± 3 


200° 


70 


Toscano et al. (1999a) 


J1730-2304 


20.5 


± 


0.4 






510" 


>50 


Toscano et al. (1999a) 


J1744-1134 


18.64 


± 


0.08 


-10.3 


± 0.5 


360P 


31 


Toscano et al. ( 1 999a) 


B1937+21 


-0.130 


± 


0.008 


-0.469 


± 0.009 


3600" 


87 


Kaspi. Taylor. & Rvba (1994) 


J 1944+0907 


12.0 


± 


0.7 


-18 


± 3 


1800" 


173 


Champion et al. (2005) 


J2124-3358 


-14 


± 


1 


-47 


± 1 


270" 


48 


Toscano et al. (1 999a) 


J2322+2057 


-17 


± 


2 


-18 


± 3 


790" 


79 


Nice & Tavlor (1995) 


Binary MSPs 


J0437-4715 


121.438 


± 


0.006 


-71.438 


± 0.007 


140P 


84 


van Straten et al. (2001) 


J06 13-0200 


2.0 


± 


0.4 


-7 


± 1 


1700" 


60 


Toscano et al. ( 1 999a) 


J0751+1807 


MA =0.35 


± 


0.03 


lip =-6 


± 2 


1150" 


22 


Nice et al. ($0?) 


J1012+5307 


2.4 


± 


0.2 


-25.2 


± 0.2 


840° 


107 


Lange et al. (20011 


J 1045-4509 


-5 


± 


2 


6 


± 1 


1940" 


119 


Toscano et al. (1 999a) 


J1455-3330 


5 


± 


6 


24 


±12 


530" 


71 


Toscano et al. ( 1 999a) 


J 1640+2224 


1.66 


± 


0.12 


-11.3 


± 0.2 


1160" 


67 


Lohmeretal. (2005) 


J1643-1224 


3 


± 


1 


-8 


± 5 


2320" 


96 


Toscano et al. ( 1 999a) 


J 1709+23 13 


-3.2 


± 


0.7 


-9.7 


± 0.9 


1390" 


57 


Lewandowski et al. (2004) 


J1713+0747 


4.917 


± 


0.004 


-3.933 


± 0.010 


1 100P 


30 


Splaver et al. (2005) 


B1855+09 


-2.94 


± 


0.04 


-5.41 


± 0.06 


910P 


17 


Kaspi. Tavlor. & Rvba (1994) 


J 1909-3744 


-9.6 


± 


0.2 


-35.6 


± 0.7 


820P 


131 


Jacobv et al. (2003') 


J1911-1114 


-6 


± 


4 


-23 


±13 


1220" 


128 


Toscano et al. (1999a) 


B 1953+29 


-1.0 


± 


0.3 


-3.7 


± 0.3 


4610" 


128 


Wolszczan et aT(2000b) 


B 1957+20 


-16.0 


± 


0.5 


-25.8 


± 0.6 


2490" 


325 


^rzournanian^t^i (1994) 


J2019+2425 


-9.41 


± 


0.12 


-20.60 


± 0.15 


1490" 


142 


Nice. Splaver. & Stairs (200 1 ) 


J205 1-0827 


1 


± 


2 


-5 


± 3 


1040" 


42 


Stappers etal. (1998) 


J2129-5721 


7 


± 


2 


-A 


± 3 


1340" 


48 


Toscano et al. ( 1 999a) 


J2229+2643 


1 


± 


4 


-17 


± 4 


1440" 


130 


Wolszczan et al. (,2000b) 


J2317+1439 


-1.7 


± 


1.5 


7.4 


± 3.1 


820" 


20 


Camilo et al. (1996) 



PDistance from parallax. 
"DM distance from NE2001. 

°Some other method used to acquire distance. The text of the cited reference should be consulted for details. 



as they oscillate in the Galaxy's potential. Our determination 
that the two velocity distributions are in fact indistinguishable 
makes this scenario unlikely. However, with identical velocity 
distributions, a difference in intrinsic luminosity distributions 
would cause the less-luminous population to be detected only 
to smaller distances and hence only to smaller scale heights. 
In fact, Bailes et al. ( 1997) find that luminosities of isolated 
and binary MSPs are different at the 99.5% confidence level, 
with the isolated MSPs being intrinsically dimmer. We have 
confirmed their results with an updated catalog; also, a simple 
examination of the median distance of the isolated population 
(510 pc) compared to the median distance of the binary popu- 
lation (1 155 pc) suggests that the isolated MSPs must be less 
luminous. 

6. CORRECTIONS TO P 

Using our upper limit on the proper motion, we can calcu- 
late the upper limit of the Shklovskii correction to the period 
derivative P. We find that 4.4 x 10" 22 or not quite 5% of the 
measured P may be due to proper motion. The acceleration to- 
ward the disk of the Galaxy is about the same size in the other 
direction, -5.0 x 10~ 22 . The acceleration in the disk makes 
a much smaller contribution, -2.2 x 10~ 23 . Combining these 
corrections changes the measured P by less than 1%. 



7. THE LOCAL INTERSTELLAR MEDIUM 

The parallax measurement of 3.3 ±0.9 mas yields a dis- 
tance of 300 ± 90 pc. This falls in between the distance es- 
timates made by the old (Taylor & Cordes 1993, hereafter 
TC93) and new (Cordes & Lazio 2002, hereafter NE2001) 
distance models: 230 and 317 respectively. Both estim ates are 
based on a DM of 4.3328 ± 0.00020cirT 3 pc JLommen et alJ 
2000). In the second galactic quadrant, between 90° and 180° 
galactic longitude, PSR J0030+045 1 is the first pulsar with a 
measured parallax, so it provides an important check of the 
NE2001 model. In addition to a parameterized model of the 
LISM, the model adds a number of clumps and voids to the 
previous TC93 model, but nothing in the direction of PSR 
J0030+045 1 . The agreement of our distance with theirs sug- 
gests no significant clumps or voids exist along this LOS. 

Nearby pulsars with known parallaxes are very useful for 
studying the LISM. Scintillation parameters measure density 
fluctuations and have been used to map out the LISM, but 
are unable to measure de nsities. The density measu rements 
must come from parallax (Bh at. Gupta. & RaolH9 98). These 
authors model the LISM explicitly as a low density bubble 
surrounded by a shell of much higher density fluctuations, but 
more pulsars with known distances are required to confirm the 
model. PSR J0030+045 1 will therefore be a marvelous tool 
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for studying the LISM. 

8. CONCLUSIONS 

We have measured the parallax of PSR J0030+0451 to be 
3.3 ±0.5 mas. We have measured its proper motion to be 
-5.74 ±0.09 mas yr _I in the plane of the ecliptic and have 
established an upper limit on its motion out of the plane at 
10 mas yr" 1 . The is one of the lowest velocities measured for 
any pulsar and is noteworthy even within the relatively low- 
velocity millisecond pulsar population. 

Combining proper motion data from this pulsar with the 
collection of existing MSP proper motion measurements, we 
find the statistical properties of the transverse velocities of iso- 
lated and binary MSPs are indistinguishable from each other. 
We do, however, find that the average z-height of the isolated 
MSPs is half that of the binary MSPs. We suggest that a lumi- 
nosity difference between the two cla s ses of objects, such a s 
that suggested by iBailes et alJ (119971) . iKramer et all (119981) . 
and Hobbs et al] (120041) . is the simplest way to account for 



both the observed difference in z-height and the similarity of 
the velocity distributions. 
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